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"PHASED ARRAY ACOUSTIC SYSTEM FOR 3D IMAGING OF MOVING PARTS" 
Description 

Field of the Invention 

The invention relates to an ultrasound, phased array imaging system and more 
particularly, to an imaging system, which uses an ultrasound transducer array with a 2D 
5 arrangement of transducer elements, having means to form 3-D motion images of moving 
parts of a body. These moving parts are typically blood flows and moving tissues like arterial 
or cardiac walls. 

The invention particularly finds applications in the field of medical ultrasound 

imaging. 

10 Background of the Invention 

Various methods and devices have been proposed for ultrasonically scanning a 
volume within a subject for three-dimensional imaging and display. Many of these 
techniques involve the scanning of a number of spatially adjacent image planes of a region of 
the body. The ultrasonic information from these associated planes can be analyzed and 

15 displayed on the basis of spatial coordinates of the data within a plane, and on the basis of the 
spatial relationship of each plane to the others. The information can be displayed in a three- 
dimensional image format such as a perspective view of the volume being imaged. A number 
of scanning techniques utilizing specially devised scanning devices have been proposed for 
acquiring these spatially related image planes. It is preferable, however, to be able to acquire 

20 information data for three-dimensional presentation without the need for special scanning 
devices or apparatus. 

It has therefore been proposed to use ultrasound phased array imaging systems 
comprising a transducer array probe with transducer elements disposed in a 2-D plane for 
forming 3-D images of a region of a body. These known systems can provide real time 3-D 

25 gray images by acquiring data in a pyramidal volume whose boundaries lean onto the 2-D 
phased array perimeter and form angles of about 30° x 30°. These 3-D images permit of 
examining in real time a body structure comprised in this pyramidal volume. In fact, these 3- 
D gray images are realized in real time with a frame rate of about 20 Hz. These 3-D gray 
images correspond to the reference volume scanned by the 2-D phased array probe where 

30 beam formation steps are performed at the level of the probe, allowing the synthesis of four 
simultaneous receive beams. Doing motion extraction in three dimensions and in real time 
presents additional problems to solve. 
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Now, conventional 2-D imaging systems for providing real time Doppler images of 
moving parts of a body are already known. Using said 2-D imaging systems, Doppler 
imaging requires analyzing typically eight successive transmit-receive signals in order to 
extract and process flow signals. With the 3-D systems previously described, real time 3-D 
5 flow imaging would necessitate an acquisition time that would be eight times longer than the 
time necessary to acquire3-D gray images of a structure, achieving a frame rate of about 2 
Hz. Another solution for acquiring real time 3-D Doppler flow images would be to acquire a 
volume eight times smaller than the pyramidal reference volume during one cardiac cycle. In 
this case, eight successive triggered cardiac cycles would be necessary to reconstruct the 

10 reference pyramidal volume. Hence, the known systems do not permit to form 3-D Doppler 
images of fluid flow or tissue motion of a body in real time because this operation would 
necessitate multiplying the acquisition time by eight in such a reference volume or would 
necessitate scanning a volume eight times smaller. In each case, the real time quality of such 
3-D images would be unsatisfactory. 

15 Because flow imaging allows to image coronary trees or vessels of organs such as 

liver and kidneys, etc., it is crucial to provide means that allow extending ultrasound 
modality to 3-D Doppler flow imaging in real time. 

Ultrasonic images are subject to image artifacts arising from a number of sources 
such as reverberation, multipath echoes, and coherent wave interference. These artifacts will 

20 manifest themselves in various ways in the images, which can be broadly described as image 
tissue. The image tissue becomes particularly troublesome when images are presented in a 
three-dimensional format, as the three-dimensional tissue can interfere with a region of the 
body and obscure said region, which the clinician is attempting to visualize. Moreover, strong 
anatomic structures like arterial walls or cardiac walls mask the weak signals generated by 

25 blood. When imaging blood flow, these structures represent the tissue. Accordingly, it would 
be desirable to provide ultrasonic image information in a format in which tissue does not 
significantly impair the images of the body region. For example, it would be desirable to 
provide ultrasonic image information in a format in which tissue information may be 
separated from flow information. 

30 It is already known to image the body using Doppler information. Doppler 

information has been used to image the body in two distinct ways. One Doppler imaging 
technique is commonly referred to as Doppler velocity imaging. As is well known, this 
technique involves the acquisition of Doppler data at different locations called sample 
volumes over the image plane of an ultrasonic image. The Doppler data is acquired over time 
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and used to estimate the Doppler phase shift or frequency at each discrete sample volume. 
The Doppler phase shift or frequency corresponds to the velocity of tissue motion or fluid 
flow within the body, with the polarity of the shift indicating direction of motion or flow. 
This information may be color coded in accordance with the magnitude of the shift or 
5 velocity and its polarity, and usually overlaid over a structural image of the tissue in the 

image plane to define the structure of the moving organs or flowing fluids. The colors in the 
image can provide an indication of the speed of blood flow and its direction in the heart and 
blood vessels, for instance. 

A second Doppler technique is known as power Doppler. This technique is 

10 unconcerned with estimations of the velocity of motion or fluid flow. Rather, it focuses 
simply on the intensity of the received signals that exhibit a Doppler shift. This Doppler 
signal intensity can be measured at each sample volume and displayed in a color variation. 
Unlike Doppler velocity imaging, power Doppler does not present the problems of 
directionality determination and low sensitivity that are characteristic of velocity imaging. 

15 Color power Doppler simply displays the Doppler signal intensity at a sample volume in a 
coded color. Like color Doppler velocity imaging, the color power Doppler display is 
conventionally displayed with a structural B mode image to define the organ or tissue 
structure in which motion is occurring. Since the value at each sample volume can be 
averaged over time or based upon a peak value, and is not subject to the constant changes of 

20 velocity and direction which are characteristic of Doppler velocity signals, the color power 
Doppler display can be presented as a stable display of motion or flow conditions in the 
body. 

Summary of the Invention 

It is an object of the invention to provide an ultrasound phased array imaging system 
25 comprising a transducer array with transducer elements disposed in a 2-D arrangement for 
forming ultrasound 3-D Doppler images in real time, such as Doppler velocity images or 
Doppler power images in real time. The system of the invention provides said ultrasound 3-D 
Doppler velocity images or Doppler power images in real time by acquiring data in a 
pyramidal volume, called reference volume, whose boundaries lean onto the 2-D phased 
30 array perimeter of the 2-D arrangement and form angles of about 30° x 30° or more. 

The systems described as prior art do not permit to form 3-D Doppler images of a 
fluid flow of a body in real time in the volume called volume of reference, because this 
operation would necessitate multiplying the acquisition time by eight with respect to the 
acquisition time for forming 3-D gray images. These systems would only permit of forming 
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3-D Doppler images of a fluid flow of a body in real time in a volume eight times smaller 
than said volume of reference, which would again result in multiplying by eight the time 
necessary to reconstruct the 3-D image of the whole reference volume. The object of the 
invention is not only to propose an ultrasound system to form 3-D gray structural images in 
5 real time, or 2-D Doppler images in real time, but to propose an ultrasound system to acquire 
3-D data and to form 3-D Doppler images of moving parts of a body in real time. 

The ultrasound viewing system of the invention comprises means to minimize the 
acquisition time duration with respect to the known systems. The number of temporal 
samples, L e. successive signals along a transmission beam, is defined as an "ensemble 

10 length". This system comprises means to reduce to 3 or 4 the number of successive signals 
out of the eight signals previously required. This operation reduces the acquisition time, 
possibly dividing the acquisition time by more than two. This system is able to provide four 
to sixteen reception beams for one emission beam, which multiplies by one to four the 
scanned volume at each transmission. Hence, The total acquisition time gain is then a factor 

15 of two to ten. The fact that few temporal samples are used is compensated by the fact that a 
larger number of spatial samples is available. A target receives an only emission beam and is 
examined using more than one received beams, for example four to sixteen received beams in 
parallel, which see this target differently, and with more than two temporal samples, for 
example three to four successive transmissions, which allows to analyze the temporal 

20 variations of the successive signals, due to the displacement of the target. 

It is also an object of the present invention to provide means for separating tissue 
information signals from fluid flow information signals. According to the invention, this 
system has means for acquiring echo signals from this reference volume and forming 
complex data. This system has means for separating tissue complex data from flow complex 

25 data. This system has further first means for tissue motion estimation and second means for 
flow estimation. The system of the invention can provide on the one hand real time 3-D 
Doppler flow velocity images or real time Doppler power images from flow estimation, and 
it can provide on the other hand real time tissue velocity images from tissue estimation. A 
technique is further presented for acquiring in real time three-dimensional ultrasonic images 

30 using 2-D phased array transducer systems. 

The invention can be applied to measure the motion rate of a moving part of the 
body, for instance the blood flow rate or tissue motion rate. 

Brief Description of the Drawings 
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The invention is described hereafter in detail in reference to the following 
diagrammatic drawings, wherein: 

FIG.l represents a general block diagram of an ultrasonic imaging system 
constructed in accordance with the principles of the present invention; 

FIG.2 represents a detailed block diagram of the ultrasonic imaging system of 

FIG.l; 

FIG.3 shows a block diagram of an ultrasound apparatus comprising the system of 

FIG.l. 

FIG.4A shows an example of transmit beam geometry; FIG.4B shows an example of 
receive beam geometry with respect to the transmit beam geometry with corresponding 
possible positions of motion estimations and FIG.4C shows other possible positions of 
motion estimations. 

Description of the preferred Embodiments 

The invention relates to an ultrasound, phased array imaging system and, more 
particularly, to an imaging system, which uses a transducer array with transducer elements in 
a 2-D arrangement. For example, the transducer elements are disposed in a 2-D arrangement 
in a 2-D plane. This ultrasound viewing system has means to form real time 3-D Doppler 
color images of moving parts of a body, for instance for imaging a fluid flow, such as a blood 
flow in a vessel or in the heart; or for imaging tissue motion such as heart walls. 

This ultrasound viewing system has means to divide by two to ten the acquisition time 
duration with respect to conventional ultrasound systems. This ultrasound viewing system 
comprises means to multiply by four to sixteen the number of beams in the reception phase. 
Hence, said system provides more than one reception beams, for example four to sixteen 
reception beams for one emission beam, which multiplies by one to four a scanned volume at 
each transmission with respect to conventional ultrasound systems. 

Besides, more than two temporal samples, for example only three or four or more 
successive signals out of eight signals are used to measure fluid flow or motion of a body 
part. This operation of using only said three or four or more successive signals again divides 
the acquisition time at least by eight to three. The fact that few temporal samples are 
available is compensated by the fact that a large number of spatial samples is available. The 
number of temporal samples is defined as "ensemble length", also called packet size, denoted 
by EL. Hence, it is an object of the invention to divide the number of temporal samples EL. 
A target receives an only emission beam and is examined using at least four and preferably 
sixteen received beams in parallel, which see this target differently, and with at least three or 
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four successive transmissions, which allows analyzing the temporal variations of the 
successive signals, due to the displacement of the target. 

The ultrasound viewing system of the invention has separation means to separate 
tissue information from flow information. Then, the ultrasound viewing system of the 
invention has further means to provide separately Doppler tissue velocity images and 
Doppler flow images or Doppler power images. A technique that is used in this separation 
means of this ultrasound imaging system is called hereafter "multiline data whitening". 

The beams Doppler signals are unknown linear combinations of the Doppler echoes 
of a moving part, such as the echoes of fluid flow, and the Doppler echoes of tissue, these 
latter echoes being temporally uncorrelated, but spatially correlated. According to the 
invention, it is possible to extract flow Doppler signals using only an ensemble length of 
three, four or more temporal samples, from spatial combination of ultrasound signals. The 
invention provides means for ultrasound "signal whitening", which is applied to 3-D flow 
imaging, allowing the separation of flow signal from tissue signal. In the example of blood 
flow estimation, the tissue is constituted for instance by vessel walls, heart walls, valves, etc.. 
A spatial correlation matrix is computed, and diagonalized allowing the separation of the 
uncorrelated Doppler components corresponding to flow and tissue signals. The Doppler 
flow signals are then processed classically with classical autocorrelation processing means to 
measure fluid velocity, for example blood velocity. The power of this signal can be used to 
provide 3-D angiograms or other cardiography images. 

Referring first to FIG.3, a block diagram of an ultrasonic diagnostic imaging system 
constructed in accordance with the principles of the present invention is shown. An ultrasonic 
probe 10 includes a 2-D phased array transducer 12, which transmits waves of ultrasonic 
energy into the body of a patient and receives ultrasonic echoes returning from structures in 
the body. In the case of ultrasonic wave transmission for Doppler interrogation of the body, it 
is the echoes returning from moving tissue, or from blood and other fluids in the body that 
are of interest. The ultrasonic probe 10 is connected to a transmitter/receiver 14, which 
alternately pulses individual elements of the transducer to shape and steer an ultrasonic beam, 
and receives, amplifies and digitizes echo signals received by the transducer elements 
following each pulse transmission. The transmitter/receiver 14 is coupled to a beamformer 
16, which controls the times of activation of specific elements of the transducer 12 by the 
transmitter/receiver 14. Circuits that are not represented, for performing transmitter/receiver 
functions and beamformer function, may be partially in the probe and partially outside the 
probe, thus forming a beamforming system 10, 12, 14, 16 in the imaging system. The timing 
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of activation enables the transducer 12 to transmit a shaped and focused ultrasound beam in a 
desired direction. The beamformer 16 also receives the digitized echo signals produced by 
the transmitter/ receiver during echo reception and appropriately delays and sums them to 
form coherent echo signals. The echo signals produced by the beamformer 16 are coupled to 
a B mode processor 61 and to the I,Q demodulator 1 8. 

According to the invention, for each transmission beam, this beamforming system 10, 
12, 14, 16 simultaneously forms in 3-D at least four to sixteen receive beams, or more, and 
respectively issues four to sixteen or more receive echo signals denoted by 4xRF to 16xRF. 

The B mode processor 61 processes the amplitude information of the echo signals, 
denoted by RF, on a spatial basis for the formation of structural images of the tissue in the 
area of the patient being scanned. According to the invention, the I,Q demodulator 18 
demodulates the four to sixteen, or more, received echo signals denoted by 4xRF to 16xRF 
into four to sixteen quadratic components, L e. complex data, denoted by 4xI,Q to 16xI,Q, or 
more, for Doppler processing. 

Referring to FIG.l, which shows a block diagram of an ultrasound system for real 
time 3-D Doppler imaging according to the invention, the four to sixteen I,Q quadratic 
complex data, or more, issued by the I,Q demodulator 18, which are denoted by 4xI,Q to 
16xI,Q, are processed by tissue/flow separation means 30 for separating 3-D tissue I,Q 
complex data from 3-D flow I,Q complex data. The separation is done using an adaptive 
spatial tissue filtering using simultaneously the four to sixteen beams acquired in 3-D. The 
tissue/flow separation processor 30 yields I,Q tissue complex data on the one hand and I,Q 
flow complex data on the other hand. The flow complex data are processed by flow 
estimation means 40 comprising autocorrelation means, which can yield flow velocity 
variance data or power data. The tissue I,Q complex data may be processed by tissue 
estimation means 43 comprising autocorrelation means, which can yield tissue velocity data. 
The tissue/flow separation processor 30 of the ultrasound imaging system of the invention 
processes the I,Q complex data in such a way that the resulting 3-D Doppler flow data are 
produced in real time. 

Referring to FIG.3, for producing flow velocity images, the flow estimation processor 
40 comprises a Doppler shift estimation processor 41, which operates in a conventional 
manner to estimate a Doppler phase or frequency shift from the I,Q flow complex data at 
each sample volume location of the image field. The Doppler shift estimation processor 41 
operates on a number of signal samples resulting from the interrogation of each sample 
volume location by an ensemble of Doppler interrogation pulses. The sample volume values 
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are applied to a flow velocity processor 51, which maps the Doppler flow values to color 
values for display. The color values are applied to a scan converter and display processor 62, 
which spatially arranges the color values in the desired image format. The color values are 
displayed as pixels on a display 60, wherein each color represents a particular velocity of 
flow in a particular direction at that pixel location. The color flow velocity information can 
be overlaid with a structural image of the interior of the body utilizing the structural 
information provided by a 3-D B mode processor 61. This 3-D compound color image can 
show both the direction and velocity of blood flow, as well as the structure of the vessels or 
organs, which contain the flowing blood. 

The Doppler system of FIG.3 can also includes a power Doppler imaging capability. 
The power Doppler components include a Doppler power estimation processor 42, which 
estimates the Doppler signal power magnitude from the I,Q signal components at each 
sample volume location using the expression (I 2 + Q 2 ) ,/2 . The Doppler power estimates at 
each location can be processed and displayed in real time or can be averaged with earlier 
acquired power estimates for each sample volume location. In a preferred embodiment, each 
sample volume location is interrogated by a number of pulses and the estimation processor 42 
utilizes the signals obtained from all interrogations in the estimations of Doppler power at the 
sample volume locations. These Doppler power estimates are mapped to display intensity or 
color values by a color power processor 52. The display values with their spatial coordinates 
are stored in an image sequence memory 64 and are also applied to the scan converter and 
display processor 62 which spatially arranges the Doppler power display values in the desired 
image format. The two dimensional Doppler power images may then be displayed on a 
display 60 or recalled from the image sequence memory 64 for three dimensional processing 
using a peak detector 63 for maximum Doppler power intensity detection. 

User operation of the system of FIG.3 is effected through various user controls 65 
which enable the user to select the type of imaging to be performed, i. e. B mode, Doppler 
color flow velocity imaging or Doppler color power imaging, and to store and retrieve 
images from the image sequence memory 64 for three dimensional display, for example. 

Referring to FIG.l, the tissue/flow separation processor 30 comprises means for 
ultrasound signal whitening. According to the invention: 

when a target receives one emission beam, then said target is examined using four 
receive beams to sixteen receive beams in parallel; 
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an ensemble length of three to four successive transmissions performed along a 
transmission beam direction allows analyzing the temporal variations of the successive 
signals due to displacement of the target; 

the receive beams of such a same multiline ensemble see the same target differently, i. 
5 e. they see the same flow and tissue Doppler components with different perspective; 

the beam signals of flow echoes and the beam signals of tissue echoes are 
substantially temporally uncorrected, but they are spatially correlated; 

the beam signals of the flow echoes and the beam signals of the tissue echoes form 
unknown linear combinations. 
10 According to the invention, a new linear combination is provided for the tissue 

information signals on the one hand and for the fluid flow information signals on the other 
hand, such that the resulting signals of flow echoes and the resulting signals of tissue echoes 
are temporally uncorrected. The technique, which is used according to the invention for 
providing the linear combinations, is called hereafter a "multiline data whitening" technique 
15 or "ultrasound signal whitening" technique. Hence the technique of the invention implies 
simple computations that can be performed in real time. 

According to the invention, it is possible to extract flow signals with only one 
ensemble length of three, four or more from spatial combination of ultrasound signals. The 
means of the invention for "ultrasound signal whitening", which is applied to 3-D color flow 
20 imaging, allows the separation of the temporally uncorrected Doppler signal components 
corresponding to flow Doppler signals and tissue Doppler signals. The correlation matrix of 
the multiline Doppler signal set is first calculated. This matrix is diagonalized. This permits 
of computing a spatial correlation diagonal matrix allowing the separation of the temporally 
uncorrected Doppler components corresponding to flow signals and tissue signals. The 
25 diagonal matrix is denoted by D. And the eigen vectors matrix is denoted by E. The number 
of temporal samples or number of successive transmissions, which can be three, four or 
more, is denoted by N. 

For each insonification depth of the volume, these definitions lead to the following 
formulations in the case of two beams, denoted by Beam 1 and Beam 2. 
30 Beam 1 X,(T) = A ltflow S fl0 w (T) + A, , tissue Stissue (T) 
Beam 2 X 2 (T) = A 2 ,fi 0 w Snow (T) + A 2 . tissue ^tissue (T) 

This formulation can be extended to a larger number N of beams. In these formulations, T= 1 
toN. 
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Since tissue and flow signals are temporally uncorrelated, two linear combinations of 
(Xi, X2) must be found comprising a linear combination Z|(T) for the flow, and a linear 
combination Z2(T) for the tissue, leading to: 

Flow Z,(T) = W, , X,(T) + W12 X 2 (T) 
5 Tissue Z 2 (T) = W 2J X,(T) + W 22 X 2 (T) 

With [Z, * Z* 2 ] =0 

[Z, * Z*,] = 1 
[Z 2 * Z* 2 ] = 1 
and C z = I 

10 Hereafter, the Flow and Tissue signals are called "Sources". The Sources are 

independent. Since neither the Sources nor the Combinations are known, the further problem 
to solve is a "Blind Source Separation" problem. However, the Sources S formed by Tissue 
and Flow Signals, are temporally uncorrelated, but spatially correlated. Hence, the Source 
Separation Resolution leads to find W verifying: 

15 Z = WX = S 

W = A" 1 

With the condition: C z = Z Z* T = I. The solution is: 

W = D" ,/2 E* T 

C X = XX* T = EDE 
20 With ZZ* T = WC X W* T = D- ,/2 E* T EDE* T ED 1/2 = I 

Where E = orthonormal eigen vectors matrix of Cx, and D = eigen values diagonal matrix of 
C x . 

These calculations are simplified in the case when two beams are considered. In this 
case, the correlation matrix is the cross-correlation of Doppler signals between beam signals 
25 for each insonification depth. A two beams system creates a 2x2 matrix: 

Cx,,. = Z X i( T )**i(T) 
Cx2, 2 =£x 2 (T)*X 2 (T) 

T=l 

Cxu=ZX 1 (T)*X 2 (T) 

T=l 

In case of normalized signals: 
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Yi-MXj 
Y 2 =X2/||X2|| 

C = E Y iC0*Y 2 CT) 



T-l 



CN -[c* ij 
x = i±Vc * c* =i± 

E ± =1/V2 



±e 



with <|> = arg(C)/2 



E = E* T = 1/V2 



Referring to FIG.4A, 4B and 4C, a very important application of this two beams 
resolution is described. This two beams resolution may be used with any order of multiline. 
FIG.4A represents one transmitted beam of direction Oz and center TB. Referring to FIG.3B, 
to said transmit beam correspond for instance sixteen receive beams. The centers of the 
1 0 receive beams are represented respectively by: 
B11,B12,B13,B14, 
B21,B22, B23,B24, 
B31,B32, B33,B34, 
B41,B42, B43,B44. 

1 5 Indeed, with sixteen receive beams acquired simultaneously, the data may be 

simultaneously processed by couple of two receive beams in order to maintain a good spatial 
resolution of the resulting 3-D data. With sixteen receive beams acquired with a square 
scheme, Bl 1, B12, B13, B14, . . . B34, B44, as illustrated by FIG.3B, one can form couples 
such as, for example: 

20 (B 1 1 , B 12) that will provide one spatial motion estimation at the position Ml , 

(B21, B31) that will provide another spatial motion estimation at the position M8, 
and so on, so that with sixteen beams, twenty four spatial motion estimations Ml to M24 may 
be achieved. Besides, four (2x2) beams forming quadruplets, such as (Bl 1,B12,B21,B22), 
may be used to estimate nine spatial motion estimations at positions Nl to N9, as illustrated 

25 byFIG.4C. 
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Twenty- four spatial motion estimations provided by associating the beams in couples, 
and nine spatial motion estimations provided by associating the beams in quadruplets give 
thirty-three possible estimated points M and N. 

FIG.2 is a block diagram that illustrates the Source Separation Algorithm of the 
5 technique of the invention. According to this technique, in a step 31 the beam correlation 
matrix Cx is calculated. Then the eigen vectors matrix E and matrix D are calculated in step 
32. The signal S is calculated in step 33. Then separation is operated in step 34. The result 
further permits estimation of Doppler power data and Doppler flow velocity data. 

The Doppler signals are further processed classically with classical autocorrelation 
10 algorithms using processor 40, as described above. 



